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Abstract 

Parton energy loss in the hot QCD medium will manifest itself not only in leading hadron spectra but also in 
reconstructed jet productions in high-energy nucleus-nucleus collisions. In this paper we report on recent theoretical 
efforts in studying full jet observables in relativistic heavy-ion collisions by discussing the modifications of jet shapes, 
inclusive jet cross section and the vector boson accompanied jet production in the presence of the QGP-induced jet 
quenching. 
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1. Introduction 

In high-energy nucleus-nucleus collisions a deconfined matter, the quark-gluon plasma (QGP) is to be formed and 
when an energetic parton produced by the initial hard scattering propagates in this hot and dense QCD medium, it 
will interact with other partons in the hot nuclear medium and be quenched substantially |[T1 [2j [3] EJ [51 . This parton 
energy loss mechanism, or jet quenching, has been proposed to be an excellent probe of the QGP and predicted 
the suppression of single hardon productions and the disappearance of away-side dihadron correlations in heavy-ion 
collisions (HIC), which later have been confirmed by the experimental measurements at RHIC and given compelling 
evidence of the existence of a new kind of matter in high-energy nuclear collisions at RHIC [6|. Not surprisingly, the 
finding of jet quenching in heavy-ion collisions has been regarded as one of the most important discoveries made at 
RHIC and initiated many theoretical explorations and very precise experimental measurements Q. 

So far the dominant experimental measurements of jet quenching is about the production of one or two hadrons 
with a large transverse momentum, which are only the leading fragments of a jet. The parton energy loss mechanism 
will not only manifest itself in leading particle productions, but more interestingly in full jet observables. Recently in 
heavy -ion community a large amount of effort has been invested in measuring reconstructed jets in high-energy nuclear 
collisions for the first time 1 8 l9l [TOl . and theories to address the novel features of full jet tomography in relativistic 
heavy -ion collisions have emerged. In this article we review our theoretical studies on jet observables of high-energy 
nucleus-nucleus collisions in perturbative Quantum Chromodynamics (pQCD) by focusing on the intra-jet energy 
flow, the inclusive jet cross section and Z°/y* -tagged jet production IfTTl [12l [131 as well as possible non-perturbative 
effects in jet productions. 
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Figure 1: Comparison of numerical results from our theoretical calculation to experimental data on differential jet shapes at y5 = 1960 GeV by 
CDF II 1 18 ] (left panel). Numerical simulation of jet shapes in vacuum , medium-induced jet shapes, and total jet shapes in medium for Pb+Pb 
with -Jsnn = 5.5 TeV at LHC (right panel). 



2. Jet shapes in HIC 

One of the most common obesevables of resolving the internal jet structure is the jet shape, which describes the 
energy distribution in a jet as lfT4ll : 

E/(£r)i®(>-#/iet) 



Here r, R are Lorentz-invariant opening angles, Rjj et - J(yi - yj e t) 2 + (<pi - 0jet) 2 gives the distance between a parton 
i and the jet in the space of rapidity y and azimuthal <p. 

In leading order we can derive the differential jet shape as ifTTIl 



dV int (r,R) V a s 2 f lz 
^ = dr = ^ 2n ~r J d zz p a^bc(z), (2) 



where 



Z = max|z, m „, — j^j if r < {R sep - l)R , 

Z = max|z mi „, R r R j if r > (R sep - l)R . 

Here 1 < R sep < 2 is introduced to take into account features of experimental cone algorithms, employed to improve 
infrared safety. In Eq. (|2|i r — (1 - z)p is related to the opening angle p between the final-state partons, P a -+b c (z) is the 
splitting kernel in DGLAP evolution equation, z is the momentum fraction of the radiated parton relative to the parent 
parton, Zmin - E"j n jEj with Ej representing the transverse momentum of a jet, and E n ™ the minimum transverse 
momentum of a parton in this jet . 

Carrying out the integration of z we get the jet shapes at LO. For instance, the jet shape for a quark gives: 



Cf(X s 2 [ ^ 1 1 Zmin 3 

2n r 



[d-z) 2 -4j) 
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Figure 2: The mean relative jet radii (r/R) in the vacuum, with complete energy loss, and in the QGP medium for ijs = 5.5 TeV central Pb+Pb 
collisions at the LHC (left panel); The ratios of total jet shape in heavy-ion collisions at LHC to the jet shape in the vacuum (right panel). 



Eq. ([3]) is collinear divergent when r — > 0. To get rid of this divergence, we need to make Sudakov resummation. 
Also we could improve the prediction power of our analysis by estimating the non-perturbative effect and include 
contributions of the initial-state radiation and corresponding Sudakov resummation [ 14 1. Taking into account all these 
contributions, we do numerical simulation and show the comparison of the theoretical results for jet shapes in p + p 
collisions with ~\/s NN = 1960 GeV at Tevatron with CDF II data in Pig. IT] It could be seen that the theoretical model 
describes CDF II data very well, especially jet shapes away from r = [1 1 J. 

In a hot and dense QCD medium, a propagating fast parton may lose a large amount of energy by medium- 
induced gluon radiation, and some lost energy carried away by radiated gluon may be recaptured by the jet and thus 
gives additional contributions to jet shapes. In our studies, we calculate the radiated gluon angular distribution in 
medium within the GLV formalism [4| and compute the medium-induced jet shapes as well as the in-medium jet 
cross sections (see Sec|3]and Sec|4]) using full numerical evaluation of the medium-induced contribution to the parton 
showers. In the model energetic inclusive jet production and tagged jet production are rare processes that follow 
binary collision scaling ~ d 2 Nu n /d 2 x with d 2 x representing the area in the transverse plane; in contrast, the medium 



is distributed according to the number of participants density 
of the QGP incorporated. 

The total jet shapes in HIC then is given bv lfTTI : 



d Npart/d X and the longitudinal Bjorken expansion 



<At< 



(1 - e) y. 



med. 



(H 



(4) 



where (1 - f q g ) ■ e gives the fraction of the energy of the parent parton falling outside of jet area with a radius R. In 
the right panel of Fig.[T]we illustrate jet shapes in vacuum \jj mc {r/R), medium-induced jet shapes ip me d.(r/R), and total 
jet shapes in HIC iff tot.(r / R) in Pb+Pb collisions with y/spfN = 5.5 TeV. It is observed that though medium-induced 
jet shapes are quite distinct from jet shapes in vacuum, the difference between jet shapes in p+p collisions and that 
in Pb+Pb is not very large. The underlining reason for this surprising result is that although medium-induced gluon 
radiation produces a broader i// me ±(r/R), this effect is offset by the fact that jets lose a finite amount of their energy, 
which will shift jet shapes in vacuum to higher Ej with steeper profile. The left panel of Fig. [2] shows the mean 
relative jet radii (r/R) in the vacuum and in the QGP medium created at the LHC for two different cone selections 
R = 0.4 and R = 0.7. We can find that the broadening of radii (r/R) of jet shapes in HIC is modest, which implies only 
a fraction of energy of jets lost in the hot medium and the QGP is rather 'gray' instead of 'black' [11]. The deviation 
of i/t t ot.(r/R) from i/f vac ,{r/R) is pronounced in the tail of jet shapes for small cone size R = 0.4, as illustrated in the 
right panel of Fig. [2] 
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3. Inclusive jet cross sections in HIC 

Parton energy loss in the QGP not only changes the energy distribution in a jet, but also alters the inclusive jet 
spectrum in high-energy nucleus-nucleus reactions. A first study of inclusive jet productions in HIC was carried out in 
the leading order (LO) ifTTI . and recently to take full advantage of jet physics in reactions with ultra-relativistic nuclei 
a next-to-leading (NLO) calculations at 0{a]) have been made fl2l . 

In hadron-hadron collisions the inclusive jet cross section at NLO can be given as fl5l : 

do ~ iet 1 C Jlr „ do ~[1 -» 2] „ lfr ± ^ 1 f Jfr _ tl dcr[2^3] a , r . 

3) . (5) 



dEjdy 



1 r dcr[2 -» 2] If do-[2 -> 3] 

= — d{E T ,y,^>} 2 — —S 2 ({E T ,y,<p}2)+— d{E T ,y,<p} } — —-S 3 ({E T ,y,( 

2! J d{E T ,y,<p} 2 3! J d{E T ,y ,^} 3 



where E Ti ,y t , (p, are the transverse energy, rapidity, and azimuthal angle of the i-th particle (i =1,2, 3), respectively, 
and <t[2 — > 2], cr[2 — > 3] stand for the partonic cross sections. In Eq. Q 5 2, S3 give the space constraints and contain 
the information of the jet-finding algorithm. At leading order, a jet is equivalent to a parton and we always have 
S2 = Zf=i — 2?=i 6(Et ( — Ej)5(yi - y). At NLO due to higher order corrections S3 has two contributions: one 
coming from the possibility that each jet still contains one parton; another resulting from the chance that a jet contains 
two partons imposed by the jet finding algorithm CD . Numerical simulation shows that the inclusive jet spectrum in 
p+p collisions with ysjyjv = 200 GeV at RHIC is described very well by a NLO code of jet production in hadronic 
collisions lfl5l . which will later be extended to study inclusive jet productions at NLO in HIC. 

In the same spirit of Eq. Q for jet shapes in HIC, we derive the related medium-modified jet cross section as 
follows fl2l : 



(r^X d^E T dy J e=0 " If. ^ '(1 - (1 - f q , g ) ■ ef d^dy 



Here T AA is the nuclear overlap function, and <r°™ the geometrical A+A cross section 1 16]. The measured cross section 
is then a probabilistic superposition of the cross sections of protojets of initially larger energy E' T = E T /(l-(l-f^ g )-e). 

Now we are ready to compute the nuclear modification factor of jet cross section B??. defined as the ratio of jet 
cross section in HIC to that in p+p rescaled by the number of binary collisions. The numerical results for the nuclear 
modification factor of inclusive jets R 1 ^ in central Au+Au with -\Jsnn = 200 GeV at RHIC are presented in the left 
panel of Fig. [3] 1 12 1, where the band of theoretical calculations shows a calculation for a ~ 20% increase in the rate 
of parton energy loss relative to our default simulation. A continuous variation of R 1 ^ with the jet radius R is found, 
which demonstrates a unique feature of jet production as compared to the single curve of Raa for pion production 
in HIC. In the calculations several cold nuclear matter (CNM) effects such as nuclear shadowing, Cronin effect and 
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Figure 4: Experimental measurements of full jet production in HIC with ^s^n = 200 GeV at RHIC |8 10]: the nuclear modification factor of jet 
production in Au+Au collisions at different jet radii (left panel); ratios of jet cross sections with R = 0.2 to that with R = 0.4 in p+p and Au+Au 
collisions (right panel). 



EMC effect have been taken into account Ifl9ll20l . We observe that for/? < 0.2 the quenching of jets approximates 
the already observed suppression in the production rate of inclusive high-p r particles. In our theoretical calculation 
CNM effects contribute significantly to the observed attenuation at large Ej. We could suppress the effects of CNM 
by taking the ratios of jet cross sections at different cone radii, which are plotted in the inserts of the left panel of 
Kg.|3] 

To see the distinction of jet productions in different systems more clearly we plot the ratio tfcr ^ ? A Q,2) / 
for p+p collisions, central Cu+Cu and central Au+Au collisions with y% = 200 GeV at RHIC in the right panel 
of Fig. [3] . It shows that the ratio decreases as the system becomes larger, and the ratio in Au+Au at RHIC is smaller 
than those in p+p and Cu+Cu collision at the same colliding energy, which implies in HIC more contributions of jet 
yields come from large angular radiation as compared to p+p collisions. 

Recently PHENIX and STAR have measured reconstructed jets in high-energy nucleus-nucleus collisions at RHIC 
and some preliminary results ||8][T0) are illustrated in Fig. [4] From the left panel of Fig. [4] one can see that measured 
R 1 ^ changes with the jet radius and at smaller radius the nuclear modification factor of jets is smaller, which agrees 
with our theoretical prediction lfTT|[T2l shown in the left panel of Fig. [3] Moreover the experimental measurements 
also demonstrate the ratio of jet yields at different radii in Au+Au is smaller than the one in p+p (see the right panel 
of Fig. [4]) which confirms the interesting feature shown the right panel of Fig.|3]with our theoretical study. We notice 
that though our theoretical predictions of inclusive jets describe very well the overall trends of the jet measurements 
at RHIC, some deviations exist if we compare the results in Fig. [3] with the data in Fig. [4] To confront the theoretical 
studies with the experimental data on a more solid base, we should include the correction from non-perturbative 
effects of jet productions ETI l22l |23| in our theoretical investigation (please see Sec|5] for more discussions), and 
precise experimental measurements of jets in HIC with large statistics will also be needed. 



4. Tagged jet productions in HIC 

A natural way to extend the study of the inclusive jet production in high-energy nucleus-nucleus reactions is to 
consider tagged jet productions in HIC, especially the jet production accompanied with a vector boson ( y or Z° ) 
in nuclear collisions. In the leading order, the vector boson and the accompanied jet are produced back-to-back in 
the center-of-mass frame and the transverse momentum of the vector boson E^ B is the same as that of the tagged jet 
E 3 T before they travel in the QGR Because the vector boson doesn't interact with the medium strongly its transverse 
momentum of the vector boson is intact. By measuring £^ B and E 3 £ in the final stage of heavy-ion collisions we may 
know exactly the energy loss of the jet by assessing E^ B - E S j. Therefore the vector boson + jet production has been 
considered as an optimal channel to pinpoint the parton energy loss in the QGP and reveal the detailed information of 
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Figure 5: Comparison of the theoretical calculation of ZP/y* tagged jet production in LO and NLO with MCFM to experimental results from 
the Fermi Lab Tevatron Collider 1 26 1 for the cross section of jets associated with Z°/y* — > p + + p~ in p + p collisions at ^ = 1 .96 TeV (left 
panel). The NLO p^-differential cross section per nucleon pair for jets tagged with Z/y* — > p + + fi~ in p+p and central Pb+Pb reactions when 
92.5 GeV < p£ l < 1 12.5 GeV (right panel). 



the hot and dense QCD medium formed in relativistic HIC. However, higher order corrections will break the momenta 
balance between tagged jet and the vector boson due to additional radiation. Here we will review our recent study of 
Z°/y*+jet production in Pb+Pb reactions at NLO and for a more detailed discussion we refer to |[T3ll24ll . 

For the baseline evaluation of tagged jet cross sections at NLO in hadron-hadron reactions we use the Monte Carlo 
for FeMtobarn (MCFM) numerical code ll25l . The confrontation of theoretical computation of tagged jet at MCFM 
NLO against the experiment measurement for p + p collisions at Tevatron in the left panel of Fig. [5] attests to the 
validity of pQCD calculations. In heavy-ion collisions the medium-induced parton energy loss will modify the tagged 
jet spectrum and the quenched jet cross section is given as follows: 

d^ A = y r dep 1 do«<* 1 Pq \ 

(PVwiPPq J e ** [1 - (1 - / 9 ,>mi„, R)W d 2 p {Z) d 2 P(jei) \ [1 - (1 - / M )e)] j ' 

with pg = P(j e t)(l - (1 - fg,g) £ )- Eq. (|7j) implies that the observed tagged jet cross section in A+A reactions is a 
probabilistic superposition of cross sections for jets of higher initial transverse energy. This excess energy is then 
redistributed outside of the jet due to strong final-state interactions. Here c/cr 4 "? / ^ 2 P(z)^ 2 P(jet) are the differential cross 
sections for away-side quark and gluon jets, respectively. 

In some experiments the Z°/y*+jet final-state channel is measured without placing restrictions on the momentum 
of the vector boson l26ll . In this case, we can integrate over P(z) in Eq. ^ and derive: 

We notice that the integrated tagged jet cross section dcr AA / d 2 pQ usually unravels the same physics as the suppression 
of the cross section in inclusive jet production discussed in Sec. [3] 

To demonstrate nuclear effects for tagged jet cross sections in HIC we define I^ A for tagged jet as: 



7jet, ff/ , d(T AA dcr pp 

l AA {R, W min ) - geo / ■ (y) 

{Taa)o\ a d PT {Z )dp T (® dp T (Z )dp T (jet) 



While at LO I AA is limited to the Z°/y* trigger momentum range, at NLO it will show rich features as we will discuss 
below. 

A numerical simulation of Zr/y* tagged jet productions in p+p and Pb+Pb with s/snn = 4 TeV for a cone size 
R = 0.2 is provided in the right panel of Fig. [5] . It is found because part of the jet energy is redistributed outside 
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Figure 6: Magnitude of the hadronization and underlying event effects used to correct the inclusive jet cross section measured by CDF with jet 
radius R = 0.7 [21 1 (left panel). Modification of pr of jets due to non-perturbative effects in p + p collisions at RHIC with Pythia8.142 [27] (right 
panel). Here for all plots the corrections are taken from hadron to parton level. 



of the jet cone the pj spectrum of tagged jet production in HIC is downshifted toward smaller transverse momenta. 
More interestingly one can observe a sharp transition of I 3 ^, from tagged jet enhancement below pr(Z) to tagged jet 
suppression above pj(zy This striking feature provides a unique prediction of jet quenching for tagged jets and will be 
experimental evidence for strong final-state interactions and parton energy loss in the QGP. Please note that while the 
same energy redistribution occurs for inclusive jet production discussed in Sec. [3] the monotonically falling spectrum 
prevents the observation of such inclusive jet enhancement. 

5. Non-perturbative effects in jet productions 

In perturbative QCD calculations of inclusive jet yield and tagged jet production, a jet is the combination of a few 
partons defined by jet algorithms and what we calculate is the parton-jet. However, in experiments what we measured 
is the hadron-jet since the final-state particles are hadrons. For jet production at high transverse momentum in a hard 
process with a large momentum transfer, the difference between parton-jets and hadron-jets may be small due to local 
parton-hadron duality (LPHD) lfT5ll2TI . Nevertheless, to make a more reliable comparison between QCD theory of 
jet productions and the jet measurements, the theoretical result of jet cross sections with pQCD should be corrected 
to hadron level, or on the contrary the experimental data should be corrected to parton level lf2TI l27l l28l . 

To include these corrections two important non-perturbative effects should be taken into account: the hadronization 
effect which results in the energy loss of a jet because when partons in a parton-jet hadronize, some hadrons may fall 
outside of the jet ( a "splashing-out" effect ); the underlying event effects which gives the energy gain of a jet due to 
interaction of beam-beam remnants not associated with the hard scattering ( a "splash-in" effect ). The good news is 
that the hadronization effect and the underlying event effect go in opposite direction and their contributions will be 
partially cancelled. 

The left panel of Fig.[6]gives the magnitude of non-perturbative effects used to correct the pj spectra with a cone 
size R = 0.7 for p + p collisions with ^Sm/ = 1960 GeV by CDF at Tevatron. We find that when pj > 200 GeV 
the fractional corrections due to non-perturbative effects are within few percent. However, this observation holds true 
only for this special case and corrections from the non-perturbative effects may vary with the colliding energy -\Js NN , 
jet transverse momenta p T and the jet radius R If23l l27l . At RHIC with colliding energy 200 GeV the measured jet 
does not have a very large px, and with small jet radii (R = 0.2, 0.4 for measurements in Fig. |4]) the contributions 
from non-perturbative effects may be considerable. In the right panel of Fig. [6] we show a numerical simulation of the 
jet pj shifting [ 27 1 to inclusive jet production due to the hadronization effect and the underlying event effects in p + p 
collisions at RHIC with Pythia8.142 l29l . It could be seen that when pj ~ 12 GeV the relative shifting Spjpr due 
to hadronization could be about 8%, and its effect will be significantly amplified by the steeply falling p T spectrum of 
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jets when we studying the jet cross section. Thus it is important to take into account corrections of non-perturbative 
effects to jet productions in HIC at RHIC and see how jet spectra will be modified by these corrections [27]. 
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